Mon. Not. R. Astron. Soc. 000, ??-H3l (20101 Printed 21 September 2011 (MN I^T^ style file v2.2) 



Asteroseismology of old open clusters with Kepler: 
direct estimate of the integrated RGB mass loss in 

. . . NGC6791 and NGC6819 

^ ! 

^ ■ A. Miglio^*, K. Brogaard^, D. Stello^ W. J. Chaplin\ F. D'Antona^ J. Montalban^ 
Oh S. Basu^ A. Bressan^'^ F. Gmndahl^ M. Pinsonneault^°, A. M. Serenelli^\ 
ly^ Y. Elsworth^ S. Hekker^^'\ T. Kallinger^^, B. Mosser^"^, P. Ventura"^, A. Bonanno^^, 

A. Noels^ V. Silva-Aguirreis, R. Szabo^^ J. Li^^, S. McCauliff^^ C. K. Middour^s, 
^ H. Kjeldsen^ 

1 I ^School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom 
Ql^ ' ^Department of Physics and Astronomy, University of Victoria, P.O. Box 3055, Victoria, BC V8W 3P6, Canada 

I ^ Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia 

• , '^INAF-Osservatorio Astronomico di Roma, via Frascati 33, Monteporzio Catone (RM), Italy 

i-S^ ■ ^ Institut d'Astrophysique et de Geophysique de I'Universite de Liege, Alice du 6 Aout, 17 B-4000 Liege, Belgium 

Q-C '^Department of Astronomy, Yale University, P.O. Box 208101, New Haven CT 06520-8101, USA 

SISSA, via Bonomea, 265, 34136 Trieste, Italy 

^ INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 

. ^Department of Physics and Astronomy, Building 1520, Aarhus University, 8000 Aarhus C, Denmark 
CZ? , 10. 



I 

O 

5h 



'Ohio State University, Dept. of Astronomy I40 W. 18th Ave. Columbus, OH 43210 USA 
■ Institute for Space Sciences (CSIC-IEEC), Facultad de Ciencies, Campus UAB, 08193 Bellaterra, Spain 

, ^■^Astronomical Institute Anton Pannekoek, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands 
^— H ' ^'"^ Department of Physics and Astronomy, University of British Colombia, 6224 Agricultural Road, Vancouver, BC V6T IZl, Canada 
^ ■ LESIA, CNRS, Universite Pierre et Marie Curie, Universite Denis Diderot, Observatoire de Paris, 92195 Meudon cedex, France 
^^INAF-Osservatorio Astrofisico di Catania, Via S.Sofia 78, 95123, Catania, Italy 
, Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748, Garching bei Miinchen, Germany 

^"^ Konkoly Observatory of the Hungarian Academy of Sciences, Konkoly Thege Miklos ut 15-17, H-1121 Budapest, Hungary 
^ 1&SETI Institute/NASA Ames Research Center, Mofjett Field, CA 94035 

Orbital Sciences Corporation/NASA Ames Research Center, Mofjett Field, CA 94035 



k> , ABSTRACT 

■ Mass loss of red giant branch (RGB) stars is still poorly determined, despite its crucial 

' role in the chemical enrichment of galaxies. Thanks to the recent detection of solar- 

like oscillations in G-K giants in open clusters with Kepler, we can now directly 
determine stellar masses for a statistically significant sample of stars in the old open 
clusters NGC6791 and NGC6819. The aim of this work is to constrain the integrated 
RGB mass loss by comparing the average mass of stars in the red clump (RC) with that 
of stars in the low-luminosity portion of the RGB (i.e. stars with L < L(RC)). Stellar 
masses were determined by combining the available seismic parameters t'max and Av 
with additional photometric constraints and with independent distance estimates. We 
measured the masses of 40 stars on the RGB and 19 in the RC of the old metal-rich 
cluster NGC6791. We find that the difference between the average mass of RGB and 
RC stars is smah, but significant (AM = 0.09±0.03 (random) ±0.04 (systematic) Mq). 
Interestingly, such a small AM does not support scenarios of an extreme mass loss for 
this metal-rich cluster. If we describe the mass-loss rate with Reimers' prescription, a 
first comparison with isochrones suggests that the observed AM is compatible with a 
mass-loss efficiency parameter in the range 0.1 ^ ?7 ^ 0.3. Less stringent constraints on 
the RGB mass-loss rate are set by the analysis of the ^ 2 Gyr-old NGC6819, largely 
due to the lower mass loss expected for this cluster, and to the lack of an independent 
and accurate distance determination. In the near future, additional constraints from 
frequencies of individual pulsation modes and spectroscopic effective temperatures, 
will allow further stringent tests of the Aj/ and fmax scaling relations, which provide 
a novel, and potentially very accurate, means of determining stellar radii and masses. 
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Figure 1 ■ Colour— maRnitu de diagram of NGC67 9 1 (lef t panel) and NGC6819 (right panel). Photometric data are taken from 
IStetson, Bruntt fc Grundahll {2003} and iHole et"aL I 1 I2OO9I) . respectively. RGB stars used in this work are marked by open squares 
and RC stars by open circles. See Sec. 13.21 for a description of the target selection. 



1 INTRODUCTION 

The mass-loss rates of evolved stars are of primary impor- 
tance for stellar and galactic evolution models, but neither 
the observations nor the theory are adequate to place reliable 
direct quantitative constraints. The mass loss from Asymp- 
totic Giant Branch (AGB) stars is reasonably well under- 
stood in terms of global pulsations lifting gas out to dis- 
tances above the p hotosphere where dust forms (e.g. IWoodI 
I1979I : iBowenlllQSa ). and the action of radiation pres sure on 
the d ust that further drives dust and gas away (e.g. iHofneJ 
l2009i . and references therein). In contrast, there is no re- 
liable theory for the mass loss of cool, dust-free red giants, 
and even the wind acceleration mechanism remains unknown 
for the Red Giant Branch (RGB). For these stars no known 
mechanisms can yet satisfactorily explain the observed wind 
characteristics (for further discussion, see lLafon fc Berruveij 
ll99ll : lHarpe"illl996l ). 

Indirect information is, however, available on the inte- 
grate d mass loss on the first giant branch (see e.g. ICatelanI 
I2OO9I . for a recent review). This is important because it 
may involve a substantial amount of mass, especially for 
the lower mass stars undergoing the helium core flash, and 
thus may affect the initial to final mass relation and the 
amount of mass recycled into the interstellar medium. In 
the case of Globular Cluster (GC) stars, mass loss along the 
RGB, and its spread, affect the distribution of stars along 
the Horizontal Branch (HB). This is because the morphol- 
ogy of the core helium burning track depends on the ratio 
between the helium core mass and the hydrogen-rich enve- 
lope mass, or, in other words, on the relative roles played 
by the core helium burning and the hydrogen shell burning 



(e.g. ISweigart fc Gross|[l976l ). The indirect evidence tells us 
that in GCs more mass is lost on the RGB than the AGB, 
since the typical stellar mass on the HB is ~0.6 Mq, with 
turn-off masses of ~0.8 Mq, and white dwarfs masses of 
-0.5-0.55 M0 (|Kalirai et al.ll2009D . 

The standard picture that the mass distribution along 
the HB matches mass loss on the RGB has been questioned 
in the last decade. The ubiquitous presence o f multiple pop- 
ulations in GCs has now beco me evident (e.g. iGratton et al.l 
I2OOII : ICarretta et"ai, 2009; Piotto et al.ll2007l ). and the pro- 
posal that the stellar helium content of the different popu- 
lations also has a role in determining the HB m orphology 
iD'Antona et al.l |2002| : iD'Antona fc Caloil |2004 ) was con- 
firmed by the pre sence of mu ltiple main sequences in some 
globular clusters (|Piotto et al., 2005 . 20(371 ) 

In order to describe the mass loss along t he RGB, 
researchers rely on empirical laws like that of iReimerj 
(1975i0), which is based on observations of Population 
I giants, and which describes mass-loss rates as a func- 
tion of stellar luminosi ty, radius, and surface gravity . 



While subsequent wor k jMullarj Il978l: iGoldberg. l1 Il979l: 
Judge fc Stencell Il99ll : ICatelanI I2OO0I : ISchroder fc Cunt j 
20071) led to slight refinements, a variant of the "Reimers' 



law" is generally used to compute stellar evolution models 
of cool stars at all ages and metallicities. 

Several observational indi cations on the RGB mass loss 
have emerged in recent years. lOrigha et all (|2002l ) detected 
the circumstellar matter around GC red giants using the 
ISOCAM camera onboard the ISO satellite, and derived 
mass-loss rates for a few RGB stars. A first e mpirical mass- 
loss law was proposed bv lOriglia et al.l (|2007l ). based on the 
observations by IRAC, onboard Spitzer, of the GC 47 Tuc. 
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They found mass loss to be episodic, depending less on 
luminosity than in Reimers' law, and occurrins predom i- 
nantly in the upper ^2 mag of the RGB (OriElia ct al. 201u). 
These findings were que stioned by [Boycr ct al. (2010.) and 
iMcDonald etall (|2011al lbh. who found no evidence for dust 
production in stars with L < 1000 L©. By comparing with 
HB models they showed that mass loss on the RGB of 47 Tuc 
is likely to be smaller than 0.24 M0. Constraints on the 
dust production and mass loss are now becoming available 
for giants belonging to other GCs jMcDonald et all |2009| . 
iBover et al.l |2009| . iMcDonald et all l201ld ." and the survey 
by Origlia et al.), and will eventually provide a complete in- 
vestigation of the differential mass loss among clusters with 
different metallicities and HB morphologies. 

Asteroseismology promises to shed new light on this 
problem by directly measuring the masses of giant stars 
in open clusters. The NASA Kepler mission, which was 
launched successfully in March 2009, includes in its field 
of view two old open clusters, NGC6791 and NGC6819. 
In these clusters sola r-like oscillations were detected in 
about 100 giants (see IStello et al.] l20ld : iBasu et all l201ll : 
iHekker et al.llioill and Fig. [l]), providing independent con- 
straints on the masses and radii of stars on the R GB, as well 
as on the distance to the clusters and their ages (|Basu et al.l 
H). 

The goal of this paper is to constrain the RGB mass 
loss by comparing the average mass of stars in the red clump 
(RC) with the low-luminosity portion of the RGB (i.e. with 
L < I/(RC)). The structure of the paper is as follows. In Sec. 
[5] we describe the procedures used to estimate the masses of 
stars from the available seismic and non-seismic constraints. 
We then address the specific case of NGC6791 in Sec. [3l 
and of NGC6819 in Sec. 2] A first comparison with model 
predictions is presented in Sec. (5] while a brief summary and 
future prospects are given in Sec. [B] 



2 ESTIMATING STELLAR MASSES 

To estimate stellar masses we use the average seismic pa- 
rameters that characterise solar-like oscillation spectra: the 
so-called average large frequency separation (Av), and the 
frequency corresponding to the maximum observed oscilla- 
tion power (t-max). The large frequency spacing is expected 
to scale as the square root of the mean density of the star: 



M/Mg 



(1) 



where Auq — 135 /iHz. The frequency of maximum 
power is appr oximatively proportional to the acoustic cut- 
off frequency llBrown et al.lll99ll:lKieldsen fc Beddindl 19951 : 



iMosser et aLll2010l : iBelkacem et al.ll201ll '). and therefore 



M/Me 



(i?/RQ)Vrcff/Teff,0 



(2) 



where i^max,© = 3100 fiRz and T^ff,© = 5777 K. 

These scaling relations are widely used to estimate 
masses and radii of red giants (see e.g . IStello et al.l l2008l : 
iKallinger et al.l2010l : [Mosser et al.l2O10l 'i. but they are based 
on simplifying assumptions that remain to be independently 
verified. Recent advances have been made on providing a 
theoretical basis for the relation between the acoustic cut-off 



frequency and z^max (|Belkacem et alJ 2011 ), and prelim inary 
investigations with stellar models |stello et al.ll2009l ) indi- 
cate that the scaling relations hol d to within y 3% on the 
main sequence and RGB (see also lWhite et al.ll201ll ). 

Depending on the observational constraints available, 
we may derive mass estimates from Equations [T] and [5] alone, 
or via their combination with other available information 
from non-seismic observations. When no information on dis- 
tance/luminosity is available, which is usually the case for 
field sta rs. Eg. [1] and [2| may be solved to derive M and R 
(see e.g. IKallinger et alllioiol : iMosser et al.ll2010l ): 
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In the case of clusters, we can use the dis- 
tance/luminosities estimated with independent methods 
(i.e. via isochrone fitting or eclipsing binaries) as an ad- 
ditional constraint. Including this information allows M to 
be estimated also from Eq. [T] or [2] alone (see Eq. [5] and [B] 
respectively), or in combination leading to a mass estimate 
with no explicit dependence on TbS (as in Eq. [7|: 
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In the following sections we use Equations [3] to [7] di- 
rectly to estimate M (and R) without adding any extra de- 
pendence on ste llar models. As illustrated in detail e.g. by 
iGai et al] l|201ll ). additional (so-called "grid-based") meth- 
ods to estimate M and R can be designed. These procedures 
are also based on Eqs.[l]and[2]but, by searching solutions for 
Al and R in grids of evolutionary tracks, have the advantage 
of reducing uncertainties on the derived mass and radius, at 
the price of some model dependence that we prefer to avoid 
in this study. 



2.1 Error estimates 

The formal uncertainties (at) on the masses (Mi) of the 
stars were used to compute a weighted average mass for 
stars belonging to the RGB and for stars in the RC: 

The uncertainties in these averages were estimated from the 
weighted scatter in the masses, i.e. 



■ = t[N 



N 



El 

\ (iv-i)Ef 



Due to the small number of stars in the samples, the correc- 
tion factor t[N — 1] was drawn from the Student's t distri- 
bution with A'^ — 1 degrees of freedom in order to assess the 
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confidence fimits correctly (see e.g. lChapHn et aTlllQQSh . We 
adopted 68% as the 1-a confidence level. 

To assess how well the formal fitting uncertainties re- 
fiected the scatter in the data we computed the ratio z — 
ajj/oM, where om is the weighted mean uncertainty esti- 
mated from the formal uncertainties on the masses, i.e. : 



Cm 



-1/2 



If mass scatter is dominated by random errors and the for- 
mal uncertainties reflect the true observational uncertain- 
ties, then z ~ 1. 

The uncertainty on AM — Mrgb — Mrc due to random 
errors was flnally computed as: 



+ ^4 



3 NGC 6791 

NGC6791 is an old open cluster, with age estimates rang- 
ing from ~7 to ~12 Gyr, depending on the choice of 
reddening and metallicity (see |Chabover. Green fc LiebertI 



ll999l:IStetson. Bruntt fc Grundahj|2003l'). M ore recent stud- 
ies (|Brogaard et al.ll201ll : iBasu et al.ll201ll ) point towards 
the lower end o f this age interval. This cluster has a so- 
lar [a/Fe] ratio l|Origlia et al.ll2006! }_ and is one of the most 
metal-rich clusters in our Galaxy (iFriel fc Janesll 19931) with 
a me l ^allicity of [Fe/H]= -|-0.3 to -1-0.5 (e.g. Brogaard et al.l 
l201ll : IOrigl ia et al.ll2006l . and references therein). 



3.1 Mass loss in NGC6791: a puzzling open 
question 

As summarised in the following paragraphs, an issue widely 
debated in the literature is whether the high metallicity of 
NGC6791 may lead to higher-than-average mass-loss rates 
along the RGB, and whether such enhanced mass loss could 
explain the very peculiar luminosity distribution of its white 
dwarfs. 

By examining the wh ite dwarf d a ta in th e deep HST 
photo metry analysed by iKing et all (|2005l '). iBedin et al.l 
l|2005[ l found a rich population of hundreds of white dwarfs. 
Surprisingly, the luminosity function of these stars showed 
two peaks, the dimmer one possibly consistent with the 
"normal" turnoff age of ~ 8 Gyr (but see below) , the brighter 
one indicating a much younger age of ~ 2.4 Gyr. White 
dwarf dating of open and globular clusters has often been 
used in the literature, and g enerally provides results consis- 
tent w ith the turnoff dating (fvon Hippelll2005l : [Hansen et al.l 



[2003). For NGC 6791. lHansenl l|2005h speculated that the pe- 
culiar luminosity function could be related to high rates of 
mass loss. If all the cluster stars formed at the same epoch 
(~ 8 Gyr ago), and if there has been significant mass loss on 
the red giant branch, many stars may have been reduced 
to a mass below the threshold to ignite the core helium 
fiash. They would evolve much more slowly than the nor- 
mal carbon-oxygen white dwarfs in the cluster, leading to 
a more luminous peak in the white dwarf lumino sity func- 
tion. This ex planation proposed b v lHanserJ ()2005l ) received 
support from iKafirai et atl l|2007l '). who derived a very low 



average mass (M~0.43±0.06 M©) for the luminous white 
dwarfs for which they obtained Keck/LRIS spectra, well be- 
low the limit for helium ignition in a degenerate core. In 
addition, the post-helium fiash, helium-core-burning stars 
in this cluster are div ided into a red gian t clump and an ex- 
tremely blue HB Ce.g. lKalirai et al.|[2007l '). The latter may be 
populated by the less massive stars evolving from the RGB 
that were still able to ignite a late helium cor e fiash (e.g. 
ICastellani fc CasteUani,1993i : lBrown et al.|[200ll '). Stars with 
slightly higher mass loss do not ignite helium and evolved di- 
rectly into helium white dwarfs, giving rise to the more lumi- 
nous peak in the white dwarf luminosity fun ction. Since the 
mass at the turnoff of NGC 6791 is ~1.1 Mr^i (jBrogaard et al.l 
a huge mass loss, possibly connected to the high 
metallicity of the cluster, must have occurred if the above 
scenario is true. 

This tentative explan at ion was disputed by 
Ivan Loon. Bover fc McDonaldl l|2008l ). who found that 
the RGB luminosity function does not show any signs of 
star depletion, and that the number of HB (clump) stars 
is consistent with most cluster stars evolving through the 
helium flash. In addition, dust - a sign of extreme mass 
loss - was not detected ar o und th e most luminous giants. 
More recently, iBedin et al.l (|2008al ) used additional data to 
examine the white dwarf sequence in NGC 6791, confirming 
a second dimmer peak in the luminosity function. Since 
the first peak is ~0 .75 mag brighter than the dimmer one, 
iBedin et al.l (|2008bh proposed that binary white dwarfs are 
responsible for the brighter peak. This explanation required 
that more than ~30% of the white dwarfs are binaries. In 
this context, the extreme HB can be attributed to binary 
evolution following the same scenario as the production 
of extreme B a nd O subdwarfs in the Galactic field (e.g. 
iHan et al.llioO^ ). 

In the new data analysed bv lBedin etall l|2008bl ). a he- 
lium white dwarf population does not seem to agree with 
the colour distribution of the cooling sequences, while a fur- 
ther problem emerged for the interpretation of the dim- 
mer peak, in that its cooling age is ~6 Gyr, compared 
with the ~8 Gyr turnoff age. This latter problem has been 
recently addressed by including in the modelling the en- 
ergy released by the sedimentatio n of ^'^Ne and of the C-0 
mixture in the white dwarf core (iDeloye fc BildstenI |2003 : 
iGarci'a-Berro et al.l l2010l : lAlthaus et al.l bOloD . which has 
the effect of lengthening the e volution of the C-O white 
dwarfs. The recent analysis by ICarcfa-Berro et al.l (1201 ll ) 
also favours a significant population of unresolved binaries 
as the most likely explanation for the observed white-dwarf 
luminosity function. 



3.2 Targets and data 

For our init i al set of targets we took the stars listed by 
ISteUo et all l|2011ah . but removed those listed as seismic 
non-members and blended stars. To avoid possible con- 
tamination by AGB stars, we then removed stars with 
L > L(RC), and subsequently separated the remaining 
stars into two groups, RGB and RC stars, by visual in- 
spection of the CMD. This classification agrees with the 
one pre sented by Stello et al . (2011a, Table 1). R ecent 
work bvlMontalban et al.l l|2010l ). lBedding et all (|201ll ). and 
[m osser et all l|201ll " has shown that the characteristics of 
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Eq. 


A/ RGB ± O"!/ 


ZRGB 


Mrc ± <y-M 


2RC 


AM ± a^-jj 




1.23 ±0.02 


1.8 


1.03 ±0.03 


1.9 


0.19 ±0.04 


Ell 


1.22 ±0.01 (0.10) 


0.7 


1.20 ±0.02 (0.10) 


0.6 


0.02 ±0.02 (0.01) 


® 


1.23 ±0.01 (0.07) 


1.1 


1.14 ±0.02 (0.06) 


1.1 


0.09 ±0.02 (< 0.01) 





1.23 ±0.02 (0.02) 


1.9 


1.07 ±0.03 (0.02) 


2.0 


0.16 ±0.03 (< 0.01) 



Table 1. NGC6791: average mass of stars in the RGB and RC estimated using different observational constraints and scaling relations 
(Equations [S] [S] (6] and[7]l. The systematic error on M due to the uncertainty on the distance is reported in brackets, z (see Sec. 12. ill was 
computed taking into account only random errors in the observables, and the correlation between BC and Tcfj. No correction to the Ai/ 
scaling was applied to clump stars. 
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1.15 ±0.03 
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0.08 ±0.04 




1.14 ±0.02 (0.09) 


0.6 


0.08 ±0.02 (0.01) 
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1.14 ±0.02 (0.06) 


1.1 


0.09 ± 0.02 (< 0.01) 





1.15 ±0.03 (0.02) 


2.0 


0.08 ± 0.03 (< 0.01) 



Table 2. As in Table[T] after applying the 2.7% suggested correction to the Af scaling for clump stars (see Sec. I3.4|| . 



the 1= \ ridge in the echelle diagram can be used as an in- 
dicator of evolutionary state of low-mass giants. A detailed 
and systematic analysis of the properties of £ = 1 modes in 
cluster stars is beyond the scope of this paper. Nonetheless, 
by visual inspection of echelle diagrams we could identify 
two misclassified stars, by looking at t he typical spacing be - 
tween ^ = 1 modes (see examples in iBedding et al.l l201lf ) . 
The first one is KIC2437103, which we now identify as a 
core-He burning star, and the second one is KIC2437589 
which has a spectrum not compatible with a RC star, and 
deserves further analysis. After this selection, we are left 
with 40 stars on the RGB and 19 stars in the RC (see Fig. 

If the mass loss is very high then the helium-burning 
stars that lost the most mass will not be in the RC, but 
rather on the corresponding flat part of the horizontal 
branch at higher Tcfi and lower L. Since our sample only 
includes RC stars, one might suspect that this introduces a 
sel ection bias. However , as shown in the membership study 
by Em ais et al.l (|201ll ). there are only 12 potential low- 
mass helium-burning stars outside the RC (disregarding the 
EHB stars, which most likely result from binary evolution). 
Among these, many are likely to be blue stragglers or non- 
members (see their Fig. 1 and Table 1) . Furthermore, work in 
progress (Brogaard et al. 2011, in prep.) shows that a theo- 
retical ZAHB that matches the RC stars is brighter than the 
potential low-mass HB stars. It al so shows that the bluest o f 
these (star 2-17, first analysed bv lPeterson fc GreenI (Il998l ) 
and often associated with the HB), is actually a blue strag- 
gler. Based on this, we believe that there are no, or at most 
a couple, helium-burning stars present with masses smaller 
than those in the RC. 

The photometric t ime-series data w ere obtained by the 
Kepler space telescope l|Koch et al.ll2O10l ) between 2009 May 
12 and 2010 March 20, known as observing quarters 1-4. We 
used the 'long caden ce' observing mode (At = 29.4 min, 
Ijenkins et al.l l2010bl ). which provided about 14,000 data 
points per star. The data reductio n from raw images t o 
the final ligh t curv es is described bv lJenkiiis et~^ (|2010al ). 
ICarcfa et all l|201ll ). and lSteUo et alj l|2011ah . A detailed de- 
scription of the asteroseismic parameters used in this work 



is given by 



iHuber et al 



SteUo et ai] (|2011bl ). who used the method of 
( 20091 ) to measure Ai^ and i^max- The aver- 
age uncertainties in Af and t'max in this sample of tar- 
gets are 1.3% and 1.7%, respectively. Other ana lyses of the 
same cl uster are give n in ou r com panion pape rs: Basu et al] 
l|201ll) . lHekker et al.l (|201ll ). and lStello etHI l|2011bD . 

T he effective t empe r atures used in this s t udy are, fol- 
lowing iBasu et al.l (I2OIII) iHekker et ah! (|201ll ). IStello et all 
(|2011a^ ■ and IStello et all (|2011bl) , b ased on the {V - K) 
colour and on the calibrations of iRamirez fc Melendej 
(2005). The value we adop ted for the reddening is E(B — 
V) = 0.16 ± 0.02 mag (|Brogaard et all I2OI1I) and for 
the m etallicity [Fe/H]=+0.3. As presented by Hekker et al.l 
the estimated random star-to-star error on Teff is 
about 50 K, which we also adopt in the following analysis. 
However, systematic errors on the Teff scale due to colour- 
temperature calibrations and reddening could inc rease the 
error budget up to ~ 110 K (|Hekker et al.|[201ll ): we dis- 
cuss in Sec. I3.4l the effect of systematic errors in Teg on our 
results. 

Finally, luminosities were deri v ed usi ng V magnitudes 
from IStetson." Bruntt fc Grundahll l|2003l ). the bolometric 
corrections by Flowed i 19961 ) ■ and the di stance modulus 
(m — M )v = 13.51 ±0.06 mag estimated bv lBrogaard et al] 
(20i3) using eclipsing binaries. 



3.3 Results 

We first applied Eq. [3] to derive masses of stars in our 
sample from fmax, A;^, and Tcff . As presented in the top 
row of Table [1] we find Mrgb = 1.23 ± 0.02 M© and 
A/rc = 1.03 ± 0.03 Mq. The mass difference we obtain is 
thus AM = 0.19 ± 0.04 Mq. The value of z = 1.8 - 1.9 
reported in the first row of Table [T] suggests that the scatter 
in the population is larger than expected from the formal 
uncertainties in the masses. This could indicate that uncer- 
tainties in the observed quantities are underestimated, pro- 
vided, of course, that there is no significant intrinsic mass 
spread within the groups of RGB and RC stars. The de- 
rived masses for RGB and RC stars are shown as a function 
of apparent V magnitude in the top row of Fig. [S] 
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Figure 2. Mass of red giants in NGC6791 as a function of visual magnitude in the RGB (left panels) and in the RC (central panels). 
Random errors on the mass estimate of each star are represented by vertical error bars. The average mass of RGB and RC stars is 
represented with full lines, and the l-cr error bars with dashed lines. Right panels show the masses of RC stars after a correction to the 
Ai/ scaling is applied (see Sec. I3.4 |l. 



To check the robustness of the mass determination, we 
also estimated masses using Eq. [S] [§] and [71 thereby includ- 
ing independent information on the distance to the cluster 
(see Sec. 13. 2p . The results are presented in Table [1] where 
we consider separately the contributions from random un- 
certainties and (in parentheses) systematic errors from the 
distance. We notice that the systematic uncertainty has little 
impact on AM (^ 0.01 Mq). The correlation between T^s 
and BC has been taken into account in the error estimates. 



If we consider only RGB stars, excellent agreement is 
found between Mrgb derived using different observables 
and scaling relations (see also Fig. [2l left panels). On the 
other hand, significant disagreement is found between Mrc 
determined with different combinations of scaling relations. 
As an additional test, and to investigate the reason for 
this discrepancy, we compared radii determined from L and 
Tcff (-Rcmd), with those derived seismically through Eq. [3] 
(J^soismo). As shown in Fig. O for RGB stars the two deter- 
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Figure 3. NGC6791: ratio between radii determined using L and T^g (-Rcmd)i and those obtained via Eq. |4] (/Jseismo)- The mass of 
each star determined via Eq.|3]is colour coded. 





Figure 4. Upper-left panel: difference between the acoustic radius of RGB and RC models of same mass and radius as a function of 
relative radius. Right panels: sound speed and temperature stratification as a function of the normalised radius in the RGB model (red 
dashed line) and RC model (blue solid line). Lower- left panel: acoustic radius as a function of mass (in solar units) for the RGB and RC 
models. 



minations of radii agree within 5% (with an average relative 
difference (i?cMD — ^?seismo)/i?CMD = —0.15%). This com- 
parison can also be carried out in terms of distance mod- 
uli, instead of radii. By coupling photometric constraints 
with i?seismo of the RGB stars, we find an average distance 
modulus (m — M)v ~ 13.51 ± 0.02 mag, where the quoted 
uncertainty represents only the standard deviation from the 
mean. This det ermination is in e xcelle nt agreement with the 
value found bv lBrogaard et all (|201ll ) using eclipsing bina- 
ries {{m-M)v = 13.51±0.06, see SecES]). The situation is 
different for RC stars, whose radii obtained from Eq. [4]are 
systematically smaller than those derived from the CMD, 
with an average relative difference of —5%. This discrepancy 
could indicate a systematic error in the scaling relations for 
RC stars. 



3.4 The Au scaling relation for RC and RGB stars 

Following the asymptotic approximation for acoustic oscil- 
lation modes, Af is directly related to the sound-travel time 
in the stellar interior. We therefore expect Av to depend on 
the stellar structure. Given that stars on the RGB have an 
internal temperature (hence sound speed) distribution sig- 
nificantly different from that of RC stars, we investigated 
whether this difference could have an impact on the mass 
determination via the Ai/ scaling relation. 

We used the code ATON 

(I Ventura. D'Antona fc Mazzitellil bOOSl ) to compute stellar 
models representative of red giants in NGC6791, i.e. 
1.2 M0 models with initial heavy-element mass fraction 
Z = 0.04, and initial helium mass fraction Y = 0.30 (see 
iMontalban et al.ll2010l . for more details on the models). The 
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evolution of the models was followed through the helium 
flash. We then considered two models with the same radius 
(within 0.1%): one on the RGB and the other in the RC. 
We c omputed radial oscilla tion frequencies with the LOSC 
code (|Scuflaire et ahllioosl ) and found that the RC model 
had a mean large frequency separation ~ 3.3% larger than 
the RGB model, despite having the same mean density. 

The average large frequency separation is known to be 
directly related to the behaviour of the sound speed in the 
stellar interior, since Ajy ~ 1/(2 Tac) where Tac is the total 
acoustic radius of the star. The latter is defined as Tac = 
tac[R), with t!ic{r) = dr' / c{r') being the acoustic radius, 
c the sound speed, and R the photospheric radius. 

The upper-left panel of Fig|4]shows the layers contribut- 
ing to the difference in the overall acoustic radius. As shown 
in the upper-right panel, the sound speed in the RC model 
is on average higher (at a given fractional radius) than that 
of the RGB model, the main reason being the different tem- 
perature profile in the two models (lower-right panel). The 
difference in acoustic radius becomes significant in regions as 
deep as m/MQ ~ 0.25, which is where the boundary of the 
helium core in the RGB model is located (lower-left panel). 
Finally, we note that while the largest contribution to the 
overall difference originates in the deep interior, near-surface 
regions {r/R > 0.9) also contribute (by 0.8%) to the total 
3.5% difference in Tac (upper-left panel of Fig. |4]). 

These findings suggest that, if we intend to use the Av 
scaling to estimate stellar mean density, a relative correction 
has to be considered when dealing with RC and RGB stars. 
This relative correction is expected to be mass-dependent 
and to be larger for low-mass stars, which have significantly 
different internal structure when ascending the RGB com- 
pared to when they are in the core-He burning phase. De- 
riving an accurate theoretical correction of the scaling, 
however, is beyond the scope of this paper, as additional 
models should be considered and detailed tests should be 
made of the internal structure of such stars. Guided by the 
representative modelling presented here, we instead adopted 
an empirical approach to correct the Av scaling. 

We determined the correction factor to the large fre- 
quency separation by minimising the average relative dif- 
ference (-RcMD — Rscismo) / RciAVi for RC stars. The correc- 
tion factor derived following this procedure is 2.7%, in good 
agreement with the correction suggested by the models. Cor- 
recting the Ai/ scaling for RC stars has a direct effect on 
Mrc obtained with different scaling relations. After we ap- 
plied the correction to the /S.v scaling, the revised average 
mass of RC stars computed with Eq. |3l [5] and [7] agreed 
within 1-cr with the value obtained with Eq. [6l which does 
not contain Av. 

We note that the uncertainties in -Rscismo and -Rcmd ini- 
ply an uncertainty of 1% in the correction factor to the 
scaling relation. This leads to an additional source of system- 
atic uncertainty in the mass estimate, which is largest (0.03 
M0) when using Eq. 3 to determine AM, and zero when 
using Eq. 6. Considering also potential differential temper- 
ature scale shifts between the RC and RGB (at the level 
of 1%), we adopted 0.04 M© as a conservative estimate of 
the additional systematic uncertainty in AM. Taking into 
account the uncertainties and scatter of the average mass 
of RC and RGB stars obtained with different scaling rela- 
tions, we derived our best estimate of the difference between 



average masses of RGB and RC stars; AM = 0.09 ± 0.03 
(random) ±0.04 (systematic). 

Finally, to check for additional sources of bias due to the 
Tcft scale, we considered effective temperatures determined 
from (V — /) colours. These T^g introduce small but system- 
atic trends in the derived quantities (e.g. mass vs. apparent 
magnitude) and worsen the agreement between masses and 
radii derived using different observational constraints. We 
find, nonetheless, that AM obtained with the different tem- 
peratures scales agree within the quoted uncertainty range. 



4 NGC6819 

NGC6819 is an open cluster of near-solar metallicity 
( Bragaglia et al.l I2OOII) estimated to be 2-2.5 Gyr old 
(|Kalirai fc Tosill2004l : lBasu et al.ll2011^ . Given that the clus- 
ter is significantly younger than NGC6791, a smaller RGB 
mass loss is expected. A s for the case o f NG C 6791, we re- 
fer tolBasu et al.l (120111). Illekker et all l|201ll ). Istello et al.] 
(|2011ah . and IStello etall (|2011bh for a description of the 
asteroseismic data available for the stars in this cluster. 

We selected targets on the RGB and in the RC fol- 
lowing the same criteria used in Sec. 13.21 and retained 
in the final list 19 RGB stars and 13 RC stars. As for 
NGC6791, we adopted the effective temperatures deter- 
mined usi ng {V — K) colours, the col our-temperature cali- 
bration of lRamfrez fc Melende j (120051). and assuming a red- 
dening E{B — V) = 0.15 mag (see iBasue t al. 2011) and 
[Fe/H]=0. 

We then performed a similar analysis as for NGC6791. 
However, since no accurate determination of the distance 
modulus was ava ilable for this cluster, we proceeded as in 
iBasu et al.l(|201ll ) and estimated the distance modulus using 
the seismically determined radii of RGB stars and photomet- 
ric constraints. We found (m — M)v = 11.80 ± 0.02 mag, 
where the quoted uncertaint y represents only the standard 
deviation from the mean fsee lBasu et al.ll201ll for a more de- 
tailed discussion of systematic uncertainties on the distance 
determination due to reddening). 

We repeated the same steps as in Sec. l3.3l to determine 
the average mass of stars on the RGB and in the RC, and 
present our results in Table |3] and Fig. [5] We notice that, 
as in NGC6791, Mrgb is independent of the scaling used, 
but significant discrepancies appear in the estimated Mrc- 
Analogously to the case described in Sec l3.4l we used 1.6 Mq 
stellar models and found that a relative correction to the Av 
scaling was appropriate, and that the radii of RC stars are 
systematically smaller than those obtained using T^s and L 
(although we note that L is not an independent quantity, as 
it is determined using Eq. |4| . 

We found that the correction factor for Au that min- 
imises the relative difference between -Rcmd and -Rscismo is 
1.9%, which is smaller than for NGC6791. Again, this is in 
qualitative agreement with the correction suggested by 1.6 
M0 models, as the differences in the sound speed profile 
between RC and RGB stars of similar luminosity decreases 
when higher masses are considered. 

By combining results from different scaling relations, 
we derive AM = -0.03 ± 0.04 Mq for NGC6819. Due to 
the lack of an independent and accurate distance measure- 
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Figure 5. As Fig.|2]but for red giants in NGC6819. 



ment, however, we consider this result less robust than the 
determination of AM in NGC6791. 



5 INFERENCES ON THE MASS-LOSS RATE 

We now compare the difference between the average mass 
of low-luminosity RGB and RC stars to theoretical predic- 
tions using a revised version of the YZVAbEI isochrones and 
stellar evolutionary sequences (|Bertelli et al.lfioOSi : Bertelli 

1 http : //stev . oapd . inaf . it/YZVAR/ 



2011, private communication). We varied the coefficient rj 
describing the efficiency of the RG B mass-loss rate, w hich 
was implemented in the models (see lBertelli et al]l 19941 ) fol- 
lowing Reimers' (1975a) empirical formulation: 

^ = 1.2710-%M-lL'■^^,ff^ (8) 

where M and L are expressed in solar units, Tcff in K, 
and t in yr. A commonly adopted value of the parame- 
ter describing the m a ss-los s efficiency is rj = 0.4, see e.g. 
iRenzini fc Fusi Peccil (|l988l ). 

We considered three scenarios: isochrones computed 
without mass loss (r; = 0), with moderate mass loss {rj = 
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Table 4. As in Table [S] after applying the 1.9% suggested cor- 
rection to the Av scaling for clump stars. 



0.35), and with an extreme value of the mass-loss efficiency 
(ji — 0.7). We notice that for these relatively young clus- 
ters, such a high mass- loss rate is not enough to suppress 
the RC and AGB phases, in contrast to old GCs (see e.g. 
iRenzini fc Fusi PeccillTgisI ). 

For NGC6791 we considered a 7-Gyr isochrone, com- 
puted with initial helium and heavy-elements mass frac- 
tions Y = 0.30 and Z = 0.04. In the case of NGC6819 
we employed a 2.1-Gyr isochrone with Y = 0.28 and Z = 
0.017. Comparisons between observed and theoretical mass- 
luminosity diagrams are shown in Fig. [S] In these plots the 
masses of individual stars were determined by Eq. |6] and 
the luminosities were estimated by combining Tag and the 
photospheric radius obtained via Eq. 3] The Reimers mass- 
loss law predicts significant mass loss only near the RGB 
tip, hence no appreciable effect is expected (AM < 0.01 
Mq) in the RGB stars we considered, which have relatively 
low luminosities. With our present measurements we have 
no evidence for the time-dependence of the mass loss, as no 
significant mass gradient is found on the RGB of NGC6791. 
Current uncertainties, however, do not allow us to constrain 
the mass loss in this portion of the RGB to better than the 
total mass loss derived. 

The comparison between observed and theoretically ex- 
pected mass difference between RGB and RC stars is pre- 
sented in Fig. [T] This first comparison shows that we can 
exclude very efficient RGB mass loss, and that the data are 
compatible with mass-loss rates described with a Reimers 
parameter smaller than ^ 0.35 and, in the better constrained 
case of NGC6791, higher than 0.1. We note that the rela- 
tion between AM and rj predicted by the models has little 
dependence on the assumed age and chemical composition 
(e.g. in the case of NGC6791 we checked the robustness of 
our main conclusion considering also isochrones of 6 and 9 
Gyr, Y = 0.28, 0.32 and Z = 0.025, 0.030). 

We finally note that in the case of NGC6819, a negative 
AM is expected from the models with little or no mass loss. 
For such a young cluster AM underestimates the integrated 
mass loss on the RGB, since the initial mass of stars in the 
RC is higher than that of stars on the RGB. In the case 
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Table 3. NGC6819: average mass of stars in the RGB and RC 
estimated using different observational constraints and scaling re- 
lations (Equations |3] El |6l and0. No correction to the Au scaling 
was applied to clump stars. 
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Figure 6. Upper panel: Masses of giants in NGC6791 derived 
from Eq. |6] as a function of luminosity. Red and green points in- 
dicate RGB and RC stars. Lines represent the 7-Gyr isochrone 
computed with different mass-loss rates. The RGB part of the 
isochrone is indicated with a dashed line, while solid lines repre- 
sent the evolution from the He flash towards the zero-age HB, and 
onwards. Lower panel: Same as the upper panel, but for stars in 
NGC6819. For this cluster, a 2.1 Gyr, solar-metallicity isochrone 
was used. 



of the 2.1-Gyr isochrone considered in Figs. [6] and [71 the 
difference in average progenitor mass between RGB and RC 
stars is ~ —0.03 Mq, hence the observed mass difference is 
compatible with no mass loss, as shown in the lower panel 
of Fig. El 



6 SUMMARY AND PROSPECTS 

We used the average seismic parameters characterising solar- 
like oscillations (fmax and Au) to estimate the masses of 
stars in the RC and on the low-luminosity RGB {L < 
L(RC)) of the open clusters NGC6791 and NGC6819. 

In NGC6791, t'max and Af were coupled with addi- 
tional constraints (photometric magnitudes, Tcff, and dis- 
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Figure 7. Upper panel: The dots connected by a dashed line 
shows the difference (AM) between the average mass of RGB 
stars (in the luminosity domain of the targets) and RC stars 
as a function of the rj parameter, resulting from a 7 Gyr model 
isochrone. The grey areas represent the l-tr region of the observed 
AM. Lower panel: Same plot as the upper panel, but for stars in 
NGC6819 and a 2.1 Gyr isochrone. 



tance modulus from eclipsing binaries) to measure stellar 
radii, and to determine the average mass of stars on the 
RGB and in the RC. When considering RGB stars, we found 
good agreement between different determinations of radius 
and mass, while in the RC a larger scatter between Mrc was 
found, together with a systematic underestimation of the 
radii obtained via Eq.|4]compared to those determined using 
photometric constrain ts and the indep e ndent distance deter- 
mination obtained bv lBrogaard et al.l (|201ll ). Supported by 
a comparison between RC and RGB models, we ascribed 
this discrepancy to the Av scaling (Eq. and argued that 
a relative correction of ~ 3% to Av between RC and RGB 
stars should be adopted for stars in this cluster. 

Taking into account the uncertainties in the observa- 
tional constraints, we found that in NGC6791 the difference 
between the average mass of RGB and RC stars is small but 
significant (AM = 0.09 ± 0.03 (random) ±0.04 (systematic) 



M0). In the context of the controversial and highly debated 
results on the mass-loss rates in NGC6791 (see Sec. l3.1|) . our 
direct estimate of AAI is incompatible with an extreme mass 
loss for this metal-rich cluster. If we describe the mass-loss 
rate with Reimers' prescription (see Eq.|8}, and assume a 7- 
Gyr isochrone, we find that the observed AM is compatible 
with a mass-loss efficiency parameter of 0.1 < 77 < 0.3. The 
low mass-loss of the RC stars, coupled with a RC mass-loss 
dispersion lower than our measurement errors, provides sup- 
port for our assumption (see Sec. l3.2|l that either no helium- 
burning stars exist with masses significantly lower than the 
measured RC stars or, if they do, they must originate from 
a very different evolution. 

Constraints on the RGB mass loss set by the analysis 
of NGC6819 are less compelling, largely because this clus- 
ter is too young to exhibit a significant mass loss (at least 
according to Reimers' formula). Moreover, no accurate and 
independent information on the distance was available for 
this cluster. 

In the future, a reliable age determination of the clus- 
ters obtained by combining all available constraints will po- 
tentially allow detailed te sts of different formulations of the 
RGB mass-loss rate (e.g. ICatelanI |200G| : ISchroder fc Cunt3 
I2OO5I , I2OO7I ), eventually leading to a better understanding 
of the physical mechanism responsible for the mass loss it- 
self. Moreover, as the duration of Kepler observations in- 
creases, it will soon be possible to study in detail the pulsa- 
tional properties of stars with luminosities higher than the 
RC. These targets could provide the link between the stars 
studied in this paper and the luminous (L > lO^L©) RGB 
and AGB giants in which a connection between mass loss 
and pulsations has been suggested by severa l studies in the 
literature (see e.g., iLebzelter fc Woodll2005l . and references 
therein). 

In this work we reached a point where possible system- 
atic biases in the Av and Vmax scaling relations may become 
the largest source of uncertainty in the determination of stel- 
lar parameters. As shown in Sec. 13.41 we have found evidence 
for systematic differences in the Ai/ scaling relation between 
He-burning and H-shell burning giants. A more systematic 
study using a larger set of stellar models is needed to ad- 
dress this issue in detail. In particular, robust model-based 
corrections to the Ai/ scaling relation would allow testing in 
detail possible systematic errors in the Ummx scaling, hence 
improving the accuracy of the mass determination. 

The frequency of maximum power will in general de- 
pend on the dynamical properties of near-surface convec- 
tive layers, and on the complex interaction between os- 
cillation modes and convection. The published relations, 
however, assume a simple scaling of i^max with the acous- 
tic cutoff frequency in the atmosphere, and simplifying as- 
sumptions about the latter. It is therefore possible that 
evolutionary-state dependent corrections may be required 
for i/max also. However, the close correspondence between 
the radius anomalies expected from Ai/ effects and those 
seen in NGC6791 suggests either that the impact of such 
corrections is modest, or that there is some effective cancella- 
tion of terms in the radius scalings. The fmax scaling should 
be investigated, but a proper formulation would involve de- 
tailed convective modelling which is outside the scope of the 
current work. One promising approach would be to search 
for residual differences in stars of known mass and radius 
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once more robust theoretical corrections to Af are applied, 
and there is work in preparation pursuing this avenue of 
approach. 

In this context, pulsating red giants with independent 
masses and/or radii (e.g. binary members and nearby gi- 
ants) would allow calibration of the and t-max scaling re- 
lations, and provide an unprecedented, and potentially very 
accurate, means of determining stellar paxameters. 

Additional observational constraints that will improve 
the precision and accuracy of our conclusions include Ke- 
pler data that will soon be available for most stars in the 
clusters. Moreover, the availability in the near future of not 
only average seismic parameters sxich as i^max and but 
also of accurate frequencies of individual oscillation modes 
in red giants, will allow more detailed tests of the scaling 
relations as well as of the internal structure of these stars. 
Finally, spectroscopic constraints on T^s will also be of great 
relevance to check the validity of the adopted T^a scale, and 
to further test the robustness of the agreement we find be- 
tween global parameters of stars determined using different 
constraints (asteroseismology, photometry, distance deter- 
mined using eclipsing binaries). 
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